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Resume 


Cette  note  presente  la  mesure  et  Panalyse  du  spectre  d’un  radiateur  Airworks.  Celui-ci 
se  rapproche  du  spectre  theorique  d’un  corps  noir  a  une  temperature  de  1073  K.  La 
mesure  nous  permet  d’e  valuer  l’irradiance  dans  la  bande  3  a  5  pm  a  partir  d’une 
mesure  d’irradiance  effectuee  avec  un  radiometre  a  large  bande. 
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1.  Introduction 


Cette  note  presente  la  mesure  et  l’analyse  du  spectre  d’un  radiateur  commercial 
Airworks  modele  C-983  (voir  figure  1).  Ces  radiateurs  produisent  d’une  fa^on  simple, 
efficace  et  economique,  un  rayonnement  infrarouge  intense.  La  mesure  du  spectre 
d’emission  et  du  niveau  d’irradiance  vise  done  a  caracteriser  le  rayonnement  emis.  A 
l’aide  du  spectre  mesure,  nous  determinons  la  fraction  de  la  puissance  totale  contenue 
dans  la  bande  3  a  5  pm  apres  propagation  sur  une  distance  de  545  m.  Nous  pouvons 
ainsi  calculer  firradiance  dans  cette  bande  a  partir  d’une  mesure  d’irradiance  effectuee 
a  une  distance  de  545  m,  avec  un  radiometre  possedant  une  bande  plus  large. 


Figure  1.  Photographie  du  radiateur. 


2.  Mesure  du  spectre 

Le  radiateur  etudie  est  de  type  PSH1000,  modele  C-983,  no  de  serie  0135.  11  est 
constitue  de  deux  filaments  chauffants  enroules  autour  d’un  cone  en  ceramique  place 
au  foyer  d’un  reflecteur  metallique  concave.  Le  diametre  du  reflecteur  est  de  40  cm. 
Lorsque  le  selecteur  de  temperature  est  place  a  la  position  HI,  les  deux  filaments  sont 
chauffes.  A  la  position  LOW,  un  seul  filament  est  chauffe.  La  temperature  est  done  la 
meme  aux  deux  positions,  mais  la  puissance  totale  emise  peut  etre  variee  d’un  facteur 
2.  Les  mesures  sont  prises  a  la  position  HI  et  lorsque  la  temperature  est  stabilisee, 
e'est-ct-dire  apres  quelques  minutes  de  fonctionnement. 
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2.1  Prise  des  mesures 


La  mesure  du  spectre  entre  0.325  et  2.5  pm  a  ete  effectuee  a  l’aide  d’un 

spectroradiometre  ASD  (Analytical  Spectral  Devices,  Inc.)  FieldSpec®  Pro  FR  qui 

collecte  la  radiation  a  l’aide  d’un  faisceau  de  fibres. 

•  L’appareil  foumit  des  valeurs  spectrales  echantillonnees  aux  1  nm. 

•  Nous  avons  mesure  le  spectre  en  pla9ant  le  faisceau  de  fibres  sans  optique  de 
collection  a  environ  4  m  du  radiateur.  Le  champ  de  vue  etait  alors  d’environ  25°. 
En  pla^ant  le  radiateur  a  cette  distance,  on  s’assure  que  les  trois  spectrometres  de 
Pappareil,  qui  possedent  des  champs  de  vue  legerement  differents,  voient  tout  le 
radiateur.  L’arriere-plan  etait  alors  negligeable  par  rapport  a  remission  du 
radiateur. 

•  Nous  avons  finalement  effectue  une  correction  logicielle  parabolique  afin  de  lisser 
les  intersections  entre  les  3  bandes  de  mesure  de  Pappareil. 

La  mesure  du  spectre  entre  2.0  et  15  pm  a  ete  effectuee  a  l’aide  d’un  spectrometre 

ABB  MR-304. 

•  L’appareil  comprend  un  detecteur  InSb  et  un  detecteur  HgCdTe  (MCT). 

•  Les  spectres  sont  donnes  par  rapport  a  a  =  I/X,  la  resolution  etant  fixee  a  4  cm*1. 

•  Le  champ  de  vue  est  de  45  mrad. 

•  Pour  limiter  le  niveau  de  signal,  Pouverture  du  spectrometre  a  et£  fix6e  a  0.8  mm. 

•  Pour  pouvoir  etalonner  Pappareil,  nous  avons  d’abord  enregistre  des  spectres  d’un 
corps  noir  (Cl  modele  SR-2-32-5A)  a  T  =  494°C,  543°C,  601°C  et  625°C.  La 
gamme  de  temperatures  a  6X6  choisie  afin  que  les  niveaux  de  signal  obtenus  soient 
du  meme  ordre  qu’avec  le  radiateur. 

•  Nous  avons  finalement  mesure  le  spectre  du  radiateur  en  pla9ant  celui-ci  a  environ 
40  cm  du  spectrometre. 


2.2  Traitement  des  donnees 

2.2.1  Etalonnage  du  spectromdtre  MR-304 

Pour  chacun  des  deux  detecteurs  (MCT  et  InSb),  nous  avons  obtenu  quatre 
mesures  d’etalonnage  Xcrt(a,  T,)  ou  T,  est  la  temperature  du  corps  noir  pour 
chacune  des  mesures.  Ces  mesures  tiennent  compte  de  la  reponse  spectrale  de 
Pappareil.  A  partir  de  la  Ioi  de  Planck  [1],  on  peut  calculer  les  valeurs 
theoriques  (relatives)  correspondantes  : 
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(1) 


ou  Ci  =  1 .4388  K/crrf 1  et  C2  est  une  constante  arbitraire.  En  supposant  que  : 

X-  (CT>  Tl  )  =  #(CT)  Ycn  fa  T,  )  +  Mp)  .  (2) 


ou  g  est  la  sensibi lite  spectrale  d’un  detecteur  avec  son  amplificateur  et  n  est 
un  niveau  de  base  independant  de  la  source,  on  peut  etalonner  V instrument  a 
partir  des  mesures  prises  aux  temperatures  T„  en  effectuant,  pour  chaque 
detecteur,  une  regression  lineaire  de  la  forme  : 

Y(a)=a(a)x(a)+b(a)  .  (3) 

Le  tableau  1  donne  un  exemple  de  la  grille  de  calcul.  Puisque  le  gain  des 
detecteurs  n’est  pas  necessairement  lineaire,  cet  etalonnage  est  valide  pour  les 
niveaux  de  signal  X(g)  de  l’ordre  de  ceux  produits  par  les  corps  noirs  utilises. 
Dans  cet  intervalle,  les  valeurs  Y(Xcn)  calculees  a  l’equation3  pour  les 
mesures  de  corps  noirs,  en  utilisant  a(a)  et  6(a)  obtenus  par  la  regression 
lineaire,  s’accordent  aux  valeurs  theoriques  Ycn  de  Pequation  1  avec  un  ecart 
relatif  inferieur  a  3%  (pour  la  plage  de  <r  tel  que  le  niveau  de  signal  des 
detecteurs  est  suffisant). 

En  appliquant  la  relation  precedente,  on  peut  done  connaitre  la  distribution 
spectrale  reelle  Yrai{(5)  correspondant  a  la  mesure  Xr(J{<5)  prise  avec  le 
radiateur. 


Tableau  1.  Exemple  du  calcul  de  a(a)  et  b(a)  pourle  ddtecteur  InSb 


T  =  494°C 

T  =  543°C 

T  -  601 °C 

T  =  625°C 

cr 

(cm-1) 

Xcn(a) 

Yen  (O) 

Xcn  (O) 

Yen  (O) 

Xcn  (O) 

Yen  (O) 

Xcn  (a) 

Yen  (O) 

a(o) 

b(o) 

3001  5 

28  732 

0  9734 

36  531 

1.3670 

51  080 

1  9464 

56.818 

2.2234 

0.0435 

-0.2573 

3003.4 

28.585 

0  9718 

36  366 

1  3650 

50.868 

1 .9440 

56.581 

2.2208 

0  0436 

-0  2558 

3005  4 

28.546 

0  9701 

36  294 

1.3629 

50  799 

1  9414 

56.494 

2.2181 

0  0436 

-0  2560 

30073 

28  387 

0  9684 

36  104 

1  3609 

50  532 

1 .9390 

56.203 

2.2155 

0  0438 

-0  2560 

3009.2 

28  171 

0.9668 

35.827 

1.3589 

50.155 

1.9365 

55.765 

2.2129 

0  0441 

-0  2568 

Note  Xcn  est  la  valeur  du  signal  donne  par  le  spectrometre  MR-304  et  Yen  est  le  signal  calculi  (unites  arbitrages) 
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2.2.2  Conversion  des  spectres  en  longueur  d’onde 


Le  spectre  mesure  avec  le  spectrometre  ASD  est  deja  exprim£  en  fonction  de 
la  longueur  d’onde,  par  contre  ceux  provenant  du  spectrometre  MR-304  sont 
exprimes  en  fonction  de  a  et  doivent  etre  convertis.  L’abscisse  des  spectres 
est  facilement  convertie  en  longueur  d’onde  puisque  X  =  1/a.  Cependant, 
l’ordonn£e  doit  egalement  etre  convertie.  En  effet,  puisque  f{X)dX  =  7(a)da 
pour  conserver  le  flux  emis  dans  des  bandes  £quivalentes,  le  spectre  en 
longueur  d’onde  est  donne,  a  partir  du  spectre  mesure  en  fonction  de  a,  par  : 


/(x)=y(o) 


da 

dX 


r(°) 

X2 


(4) 


2.2.3  Mise  a  I’echelle 

Puisque  les  spectres  sont  mesures  selon  une  echelle  relative,  il  est  n£cessaire 
d’effectuer  une  mise  a  l’echelle  pour  faire  corresponds  les  mesures 
effectuees  avec  les  deux  spectrometres.  On  multiplie  les  spectres^^)  mesures 
avec  chaque  appareil  par  une  constante  afin  que  les  courbes  se  rejoignent  aux 
intersections  et  pour  normaliser  le  maximum  du  spectre  a  une  valeur  unitaire. 
Pour  le  spectrometre  MR-304,  une  legere  correction  (~2%)  est  aussi 
necessaire  entre  les  detecteurs  MCT  et  InSb. 

2.2.4  Coupe  des  courbes 

On  coupe  ensuite  les  courbes  pour  conserver  les  points  dans  la  bande  ou  un 
detecteur  donne  est  le  plus  precis  (reponse  suffisante).  Les  sections 
conservees  sont  les  suivantes  :  avec  ASD,  de  0.35  a  2.4  \xm;  avec  InSb,  de  2.4 
a  4.9  pm  et  avec  MCT,  de  4.9  a  14  pm.  Les  vecteurs  X  etf{X)  correspondant  a 
ces  trois  bandes  sont  ensuite  mis  bout  a  bout. 

2.2.5  R6-6chantillonnage 

Afin  de  pouvoir  manipuler  plus  facilement  le  spectre,  on  re-echantillonne 
finalement  le  spectre  J{X)  a  chaque  1  nm  en  interpolant  les  valeurs  mesurees 
(qui  sont  au  depart  a  chaque  1  nm  pour  X  <  2.4  pm  puis  a  chaque  2  cm'1  pour 
X  >  2.4  pm). 

2.3  Spectres  obtenus 

A  partir  du  traitement  decrit  precedemment,  on  obtient  le  spectre  d’emissiony(X)  du 
radiateur  presente  a  la  figure  2.  Les  zones  plus  bruyantes  sont  dues  aux  bandes 
d’absorption  atmospherique  qui  reduisent  le  niveau  de  signal.  Selon  la  longueur 
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d’onde  du  maximum,  la  temperature  serait  d’environ  1073  K.  En  effet,  selon  la  loi  de 
Wien  [1]  : 


Xnux=  2897.8  T-'  [pm],  (5) 

ou  la  temperature  T  est  en  Kelvin.  La  courbe  correspondant  a  un  corps  noir  a  cette 
temperature  est  egalement  presentee  sur  la  figure  2.  Comme  pour  le  spectre  du 
radiateur,  celle-ci  est  normalisee  a  1  a  son  maximum. 


Figure  2.  Spectre  demission  normalise  du  radiateur,  ainsi  que  celui  d’un  corps  noird  T  =  1073  K. 


3.  Mesure  de  I’irradiance  a  545  m 

Pour  mesurer  l’irradiance  a  une  distance  de  545  m  du  radiateur,  nous  avons  utilise  un 
radiometre  constitue  d’une  lentille  collectrice  et  d’un  detecteur  pyroelectrique  place  au 
foyer  de  celle-ci.  Un  hacheur  place  devant  le  detecteur  a  permis  de  moduler  la 
radiation  incidente  a  une  frequence  de  300  Hz.  Une  lunette  de  visee  a  aussi  ete 
installee  a  cote  de  l’ensemble  afin  d’en  permettre  l’alignement  angulaire.  L’acquisition 
a  ete  effectuee  a  l’aide  d’un  oscilloscope  numerique.  Voici  les  parametres  importants 
de  ce  systeme  de  mesure. 

•  Lentille  collectrice  (CaF2)  : 

■  Longueur  focale  ~  250  mm;  et 

■  Diametre  ~  44  mm  (degagement  du  support). 

•  Detecteur  pyroelectrique  Molectron  P4-41  #162  : 
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■  Surface  de  detection  de  1  mm  x  1  mm;  et 

■  Le  champ  de  vue  correspondant  est  de  4  mrad. 

•  Oscilloscope  LeCroy  Waverunner  LT344L  : 

■  Base  de  temps  de  1000  points,  2  ms/div; 

■  Couplage  de  1  MQ  AC,  BWL  25  MHz; 

■  On  applique  la  fonction  “Enhanced  Resolution”  (Eres+3  bit),  ce 
qui  augmente  la  resolution  de  f  oscilloscope  a  8  +3  bit,  en  plus  de 
produire  un  filtrage  passe-bas  (largeur  de  bande  3  dB  =  400  Hz); 
et 

■  On  applique  finalement  une  moyenne  sur  50  balayages. 

Le  systeme  de  detection  a  et£  prealablement  etalonne  avec  un  laser  emettant  dans 
Tinfrarouge  moyen  (systeme  MISLED  [2],  longueur  d’onde  de  3.726  pm,  impulsions 
de  1 2  ns  avec  un  taux  de  repetition  de  20  kHz).  Notons  que  la  base  de  temps  de 
f  oscilloscope  etait  alors  differente,  mais  cela  ne  devrait  pas  affecter  la  reponse 
puisque  f  amplitude  du  signal  crete  a  crete  etait  peu  affectee  par  le  niveau  de  filtrage 
(Eres)  selectionne.  L’echelle,  puissance  moyenne  incidente  par  rapport  k  la  tension 
crete  a  crete,  est  de  2.2  mW/V  (±15%)  pour  une  frequence  de  modulation  de  300  Hz. 
La  lecture  de  puissance  etant  prise  devant  le  radiometre,  la  reponse  mesuree  tient 
compte  de  la  transmission  de  la  lentille  collectrice  a  3.726  pm. 

Pour  les  mesures  d’irradiance  (efifectuees  le  2004-10-25)  nous  avons  enregistre  une 
trace  avec  le  radiateur  en  operation  ainsi  qu’une  autre,  farriere-plan,  sans  le  radiateur 
(voir  figure  3).  Puisque  les  lames  du  hacheur  sont  plus  chaudes  que  farriere-plan,  les 
traces  mesurees  avec  et  sans  radiateur  ne  sont  pas  en  phase.  On  ne  peut  done  pas 
simplement  soustraire  les  valeurs  Crete  a  Crete.  On  soustrait  done  les  traces  puis  on 
mesure  la  valeur  crete  a  crete  sur  la  difference.  La  puissance  ainsi  mesuree  est  de 
1.3  pW  (soit  0.61  mV  crete  a  crete),  ce  qui  correspond  a  une  irradiance  de  88  nW/cm  . 
En  pla£ant  un  filtre  passe-bande  (le  spectre  de  transmission,  -3.33  a  4.54  pm,  sera 
presente  un  peu  plus  bas)  devant  la  lentille  collectrice,  nous  avons  obtenu  environ  25% 
de  cette  valeur.  Le  rapport  signal/bruit  pour  cette  demiere  mesure  est  cependant 
beaucoup  plus  faible  puisque  le  niveau  de  signal  du  radiateur  est  alors  comparable  au 
niveau  de  signal  de  farriere-plan. 
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Figure  3.  Signal  traits  obtenu  lors  de  la  me  sure  dirradiance  du  radiateur. 


4.  Analyse 

En  comparant  a  la  figure  2  la  distribution  de  puissance  mesuree  a  celle  d’un  corps  noir, 
on  constate  que  l’emissivite  du  radiateur  semble  plus  elevee  pour  les  longueurs  d’onde 
elevees.  Cela  s’explique  probablement  par  le  fait  que  la  tole  du  reflecteur  concave 
devient  chaude  et  emet  comme  un  corps  noir  a  une  temperature  de  quelques  centaines 
de  degres  Celsius.  Si  tel  est  le  cas,  la  forme  du  spectre  devrait  dependre  de  la  distance 
de  propagation  puisque  le  rayonnement  emis  par  la  tole  chaude  serait  moins 
directionnel  que  la  reflexion  sur  le  reflecteur  concave.  Observe  a  grande  distance,  le 
spectre  du  radiateur  devrait  alors  tendre  vers  celui  d’un  corps  noir  a  1073  K. 

Pour  calculer  la  forme  du  spectre  re$u  par  notre  radiometre  situe  a  une  distance  de 
545  m,  il  faut  tenir  compte  de  la  transmission  atmospherique,  de  la  transmission  de  la 
lentille  de  CaF2  et  de  la  reponse  spectrale  du  detecteur  pyroelectrique  utilise.  Pour  la 
mesure  effectuee  avec  un  filtre  passe-bande,  on  doit  aussi  tenir  compte  du  spectre  de 
transmission  de  ce  dernier.  La  figure  4  presente  les  fonctions  de  reponse  et  de 
transmission  utilisees  pour  le  calcul.  La  courbe  de  transmission  atmospherique  a  ete 
calculee  a  l’aide  de  l’outil  logiciel  MODTRAN  pour  le  modele  "mid-latitude  summer" 
tandis  que  les  courbes  de  reponse  du  detecteur  et  de  transmission  (exteme)  de  la 
lentille  sont  fixees  approximativement  a  partir  de  [3]  et  de  [4,  p.  7-16]. 
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Figure  4.  Fonctions  de  transfert  utihs6es  pour  calculer  le  spectre  tel  que  ddtectd  £  545  m  parle 

radiom&tre 


Afin  de  montrer  l’effet  de  la  propagation  atmospherique  seule,  la  figure  5  presente 
dans  un  premier  temps  le  spectre  du  radiateur  mesure  de  meme  que  celui  calcule  pour 
une  distance  de  propagation  L  de  545  m.  Pour  L  =  545  m,  il  s’agit  done  du  spectre 
incident  sur  le  radiometre.  Notons  qu’on  n’y  tient  pas  compte  de  la  baisse  de  puissance 
en  ~1/Z,2  mais  seulement  de  la  transmission  atmospherique.  Comme  nous  Pavons 
mentionne  plus  tot,  la  baisse  de  puissance  due  a  Peloignement  de  la  source  pourrait 
l^gerement  dependre  de  la  longueur  d’onde  si  le  spectre  mesur£  devait  provenir  en 
partie  de  la  tole  du  reflecteur.  La  figure  6  presente  finalement  le  spectre  tel  que  re9u 
par  le  radiometre  (avec  et  sans  flltre  passe-bande).  II  s’agit  du  spectre  calcule  a  545  m 
de  la  figure  5,  multiplie  par  les  fonctions  de  transfert  presentees  a  la  figure  4  pour  la 
lentille,  le  detecteur  et  le  flltre  passe-bande. 

A  partir  de  ces  spectres,  on  peut  evaluer  la  fraction  de  la  puissance  totale  mesuree  par 
le  radiometre  qui  est  dans  la  bande  3  a  5  pm.  Pour  le  spectre  emis  par  le  radiateur  et 
detecte  a  L  =  545  m  par  le  radiometre,  on  obtient  environ  50%  de  la  puissance  totale 
entre  3  et  5  pm.  En  effectuant  le  meme  calcul  pour  un  corps  noir  a  1073  K,  le 
pourcentage  est  legerement  plus  elev£  puisqu’on  obtient  alors  53%  de  la  puissance 
totale  detectee  dans  la  bande  3  a  5  pm.  Puisque  la  difference  est  faible,  on  considerera 
que  la  valeur  de  50%  est  la  bonne. 

Ainsi,  d’apres  les  mesures  d’irradiance  presentees  a  la  section  precedente,  ou  nous 
avons  mesure  une  irradiance  totale  de  88  nW/cm  avec  le  radiometre,  on  peut  evaluer 
que  f  irradiance  incidente  sur  le  radiometre,  entre  3  et  5  pm,  etait  d’environ 
44  nW/cm2  (soit  50%  de  firradiance  mesurde). 
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Figure  5.  Spectre  normalise  6m is  par  le  radiateur  et  spectre  calcu!6  apr6s  propagation  sur  545  m 


Figure  6.  Spectre  d6tect6  6  545  m  par  le  radiomdtre  avec  et  sans  filtre  passe-bande 


Lors  des  mesures  d’irradiance,  nous  avions  mesure  environ  25%  de  la  puissance  sans 
filtre  en  pla^ant  un  filtre  passe-bande  (3.33  a  4.54  pm)  devant  le  radiometre.  Si  on 
considere  le  spectre  detect^  multiplie  par  le  spectre  de  transmission  du  filtre  utilise 
pour  cette  mesure  (voir  figure  6),  on  obtient  egalement  un  rapport  de  25%,  ce  qui  est 
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en  parfait  accord  avec  la  mesure.  Notons  toutefois  que  la  transmission  du  filtre  semble 
non  nulle  pour  les  longueurs  d’onde  superieures  a  8  pm  (voir  figure  4).  Nous  avons 
neglige  cette  contribution  pour  le  spectre  calcule  a  la  figure  6  puisque  nous  ne 
connaissons  pas  la  forme  exacte  du  spectre  de  transmission  du  filtre.  Etant  donne  que 
1’emission  du  radiateur  ainsi  que  la  transmission  de  la  lentille  collectrice  sont  plus 
faibles  dans  cette  bande,  cela  produirait  toutefois  une  contribution  peu  importante  sur 
la  valeur  calculee  avec  filtre. 

Finalement,  connaissant  la  distribution  spectrale  a  545  m  (figure  5),  ainsi  que 
l’irradiance  dans  la  bande  3  a  5  pm  (44  nW/cm2),  on  peut  etalonner  le  spectre  en 
irradiance.  Celui-ci  est  presente  a  la  figure  7.  L’irradiance  E  (en  nW/cm2/pm)  y  est 
evaluee  a  partir  de  : 


E(X,L  =  545)  = 


44  /,nc  (A.,  L  =  545) 
I/**  (^>  L  =  545)d  X 


(6) 


ou 


finc (X, L  —  545)  =  Tatm{L  —  545)x  /(?t)  est  le  spectre  non  etalonne  de  la  figure  5. 


50 


X  [Mm] 


L  =  545  m 


8  10  12 


Figure  7.  Distribution  spectrale  d’mradiance  £  545  m  du  radiateur 


5.  Conclusion 


Nous  avons  mesure  le  spectre  d’emission  d’un  radiateur  Airworks  ainsi  que 
1’ irradiance  a  545  m.  Le  spectre  se  rapproche  de  la  distribution  theorique  pour  un  corps 
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noir  a  une  temperature  de  1073  K.  Un  tel  radiateur  peut  done  etre  utilise  comme  source 
chaude  et  puissante.  Nous  avons  calcule  qu’a  545  m,  la  puissance  dans  la  bande  3  a 
5  pm  represente  50%  de  la  puissance  totale  mesuree  a  l’aide  de  notre  radiometre. 
Ainsi,  nous  avons  pu  evaluer  la  distribution  spectrale  d’irradiance  a  cette  distance,  qui 
est  representee  a  la  figure  7. 
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